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Abstract. One of important questions concerning particles which compose the 
Dark Matter (DM) is their average speed. We consider the model of relativistic 
weakly interacting massive particles and try to impose an upper bound on their 
actual and past warmness through the analysis of density perturbations and com- 
parison with the LSS data. It is assumed that the DM can be described by the 
recently invented model of reduced relativistic gas (RRG) . The equation of state of 
the RRG model is closely reproducing the one of the Maxwell distribution, while 
being much simpler. This advantage of the RRG model makes our analysis very ef- 
ficient. As a result we arrive at the rigid and model- independent bound for the DM 
warmness without using the standard (much more sophisticated) approach based 
on the Einstein-Boltzmann system of equations. 
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1 Introduction 

The Sun is shining bright in Brazil. However, independently on geography, cosmologists say the 
Universe is dominated by a darkness. Namely, the energy balance of the present-day Universe 
shows that the relative energy densities of the Dark Energy and Dark Matter (DM) are close 
to Q\ = 0.7 and ^^5)Af ~ 0.25, respectively, while the visible (more precise, baryonic) matter is 
represented by a modest less than 5% of the total energy density [1]. 

Despite the existing variety of the models for Dark Energy, the Cosmological Constant 
(CC) A is the most natural candidate. The presence of a A-term is dictated by the requirement 
of consistency of quantum field theory in curved space. At the same time the enormous fine- 
tuning which is necessary for adjusting the value of A to the astronomical observations creates 
a longstanding CC problem (see, e.g., [21 [3] for discussion). However, in this paper we will 
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concentrate on the second dark component, which is equahy mysterious. The main candidate 
to be DM is the gas of weakly interactive massive particles (WIMPs) which could be part 
of a multiplet composition of some extension of the Standard Model of elementary particle 
physics. For example, those can be superpartners of observable particles in MSSM or in some 
supergavity model. One can find an extensive discussion of the DM issue in the books jUEKUl?] 
or in the recent reviews [HJ |9l |10] . 

In simplest terms one can describe the DM problem as follows. The astronomical obser- 
vations show that the stars and interstellar gas clouds in the spiral galaxies have the rotation 
curves different from the ones produced by gravitational field of the visible matter. The typical 
spiral galaxy has an almost flat structure, while the gravitational field is apparently produced 
by some almost spherical distribution of mass, total amount of it should be a few times greater 
than the one of the visible part. The hidden mass presumably forms a halo and is called DM. 
The main question is from what the DM is made. Obviously, the constituents of the DM should 
have properties distinct from the ones of the baryonic matter, for otherwise the two kinds of 
matter would be distributed in the same way. Furthermore, in the cosmological setting, DM is 
necessary for the cosmic structure formation. 

One can distinguish the three main kinds of DM. The first one is cold DM, e.g., formed 
by WIMPs. Another one is hot DM, which can be represented, e.g., by massive neutrinos 
(with a mass of some eV). Hot dark matter leads to the so-called up-bottom scenario, where 
structures of clusters of galaxies are formed first, while cold dark matter implies the bottom-up 
scenarios, forming first small objects of scales smaller than a galaxy. Even if the cold dark 
matter scenario seems more favored, each scenario has its own problems, with suppression (hot 
DM) or excess (cold DM) of power at small scales [11]. The intermediate scenario of warm 
DM has been invoked to solve this problem [12^ [T3\ [T^ [T6\ 117]. Warm DM may be composed, 
e.g., by relatively heavy sterile neutrino with the keV-scale mass (see, e.g., [151 [161 [HI [19] ) or 
by some other particles such as light gravitinos [20]. The structure formation in these models 
has been explored using fluid description [l7| (see also [16] and references therein) and also 
A^-body simulations methods [211 [9]. 

The concepts of hot and warm DM imply that the DM constituents are relativistic in the 
early Universe and then cool down when the Universe expands. Qualitatively similar evolution 
takes place also for the baryonic matter, however the last emits radiation and therefore cools 
down faster. Due to the growing amount of the available experimental data, the requirements 
for the cosmological models are becoming stronger and in particular it would be desirable to 
have more precise description of the expansion of the Universe. In both cases of baryonic and 
DM one needs to have a model which could take into account, in a natural way, the continuous 
evolution of the equation of state of the matter content. In other words, this should be a model 
with the radiation-like behavior in the early Universe, becoming more like a dust or like a DM 
in the later epoch. 

In the present paper we start the detailed analysis of the new, relatively simple model of 
reduced relativistic gas (RRG). This model is reproducing the equation of state of the Maxwell 
distribution with high precision [22]. The short technical introduction to the RRG model is 
postponed for the next section, but some general discussion is in order here. The RRG model 
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enables one to consider the matter content which is hot in the early Universe and continuously 
cools down when the Universe expands. Let us note that the same purpose can be achieved by 
taking, e.g., relativistic distributions for the ideal gas case, such as the Maxwell, Fermi-Dirac or 
Bose-Einstein ones. However, in all these cases the relation between the energy density and the 
pressure have rather complicated form. In particular, for the most simple Maxwell distribution 
this relation has the form of the ratio of the two modified Bessel functions [23j (see also, 
e.g., [Mj)- Needless to say that such expressions are difficult to deal with in the cosmological 
framework. It is important to remember that the equation of state is the unique relevant 
element of the matter content, as far as we are interested only in the behavior of the conformal 
factor of the metric (zero-order cosmology). In this respect the RRG model [22] provides a real 
advantage. Due to its simplicity, one can easily integrate the Friedmann equation and develop 
the zero-order cosmological model in the economic and analytic form. From the other side, as 
we have already mentioned above, the RRG model provides the equation of state which is very 
close to the one corresponding to the Maxwell distribution. Therefore the properties of the 
RRG-based cosmological model can be safely attributed to the model of the Universe filled by 
an ideal gas of massive relativistic particles. 

In the framework of the RRG model one arrives at the single-fluid cosmology which in- 
terpolates, in a natural way, between the radiation-like and dust-like regimes. Therefore, at 
the level of zero-order cosmology the RRG model is proved to be a useful tool. Even more 
important, one can generalize the RRG model by implementing the possibility of the energy 
exchange between RRG gas and other fluids by introducing viscosity (see, e.g., [25] and ref- 
erences therein). In this way we may expect to obtain the most precise zero-order cosmology 
which can be hopefully given by an analytical solution and be very useful. In future, we may 
have a chance to see the impact of different interaction types of the DM particles, e.g., on the 
cosmic perturbations spectrum. 

However, before the RRG model can be considered as a tool for creating the precise cosmo- 
logical model, it has to be submitted to another important test. The present-day cosmological 
model does not deserve confidence if it produces good results only at the zero-order level. 
So, the next step in the development of the RRG model is to check whether it can produce 
acceptable results for the cosmic perturbations. If the RRG model can pass this test, it is 
worthwhile to use it as a building block for constructing the realistic cosmological models in 
a way described above. For this reason, in the present paper we address the issue of density 
perturbations spectrum in the simple RRG model [22]. The output of the analysis of the 
cosmic perturbations in the RRG model can be compared to existing detailed description of 
the perturbations in the warm and hot DM models (see, e.g., \12\ [26] and also [HJ [9l [5l [6] 
for the review). Indeed, our purpose is not to compete with the standard results, which are 
based on the numerical solution of the much more detailed description in the framework of the 
complicated Einstein-Boltzmann system. Instead of this, we want to see whether the results 
derived within the RRG model are compatible with the standard ones and, in this way, to try 
our model. 

Let us remember that the RRG model is reproducing the ideal relativistic gas, which is 
isotropic and, therefore, can not provide full information on the motion of the DM particles. 
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In the framework of such simphfied approach we have a restricted choice of relevant physical 
observables, the most obvious is an average speed of the DM particles. Hence, our immediate 
purpose here is to establish an upper bound for the velocities of the DM constituents, both 
in the present and earlier epochs of the Universe. The upper bound for the DM velocities 
comes from the fit with the LSS data [27] • In this way, we can test to which extent the dark 
matter can be hot or at least warm. It is well known that the standard way to impose the 
bound on the warmness of the DM particles is through the analysis of cosmic perturbations 
in the Einstein-Boltzmann coupled system [M]- As we shall see in what follows, the use of 
the relatively simple RRG model enables one to achieve similar restrictions in a much more 
economic way. Thus, using this approach, we circumvent the technical difficulties related to 
the analysis of the Einstein-Boltzmann system without losing the essential features. 

The paper is organized as follows. In the next section we present a very brief introduction 
to the RRG model. The reader can find further details in [22]. In section 3 the equations for 
density perturbations and their numerical analysis are considered and in section 4 we present 
some discussions and draw our conclusions. 



2 Reduced model for relativistic gas 



The equation of state for the ideal relativistic gas of identical massive particles has been derived 
in [23]. This equation involves a ratio of two modified Bessel functions and is rather difficult 
to apply for the cosmological purposes. One can simplify things considerably if assuming 
that, instead of the Maxwell law, all particles have equal kinetic energies. An elementary 
consideration leads to the following relation between pressure P and energy density p = ne: 



P 
3 



1 



(1) 



where e = mc^/\/l — n is a number of particles per unit of volume and /3 = v/c. One can 
introduce the new notation for the density of the rest energy 



Pd 



nmc 



(2) 



where n is the number of particles for a unit of 3d volume. This density depends on the scale 
factor in the usual way 



Pd 



~^ Pd ■ 



where is the rest energy density at present, when a = oq. Using this quantity, one can 
rewrite Eq. ([T]) in the form 



1 _ lJL 

1 

p 



(3) 



Both forms ([T|) and ([3|) will be useful for us in what follows. An important observation is that 
the expression ([3]) reproduce the Maxwell-based equation of state with very good precision. 
According to the plot obtained in [22] , the maximal difference between the equations of state 
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p = p{P) in the two cases is just 2.5% of the absolute value of p and, moreover, this discrepancy 
goes to zero pretty fast in the UV, when p ^ pd- 

Let us emphasize that the difference between the equation of state which follows from 
relativistic Maxwell distribution and the equation of state in our simplified model is so small 
that it can be seen as negligible, when we use this equation of state, e.g., in the Friedmann 
equation. Therefore, the cosmological model which we are going to develop on this background, 
will be based on the following two assumptions: 1) that the massive particles (e.g. the ones 
of DM or baryonic matter) go from one equilibrium state to another in a sufficiently smooth 
way, such that the fluid composed by these particles can be described by the equation of state 
instead of the Boltzmann equation (if compared to the standard approach [26]). 2) That the 
interaction between these particles is negligible. Indeed, the main advantage of our model is 
that it enables one to introduce interactions between the particles in a very elegant way. We 
shall treat this issue in a separate paper and now concentrate on the ideal gas case. 

Since the RRG model is really close to the Maxwell distribution, in what follows we shall 
refer to the velocity of the particles in the RRG as to "average speed" . This term will help us 
to keep in mind that the results of our calculations provide the reliable information not only 
about the proper RRG model, but also about the Maxwell-distributed relativistic gas. 

Using the conservation law 

dp 3(ia 



p + P a 

and the equation of state ([3]), one can easily arrive at the RRG density scaling rule 



(4) 



p = p{z) = -^^{l + zf^ll + h^{l + z)\ (5) 

where il^j = il2)Af + ^ total relative present-day matter energy density, p^^. is the present 

day critical density and z = —l + ao/ais the red-shift parameter. The dimensionless parameter 
h shows whether the velocity of the RRG particles is large or small or, in other words, whether 
the matter is "cold", or "warm ", or "hot". In order to better understand the physical sense 
of this parameter, let us express it in two different (albeit equivalent) forms as follows: 



h=f^ = ^^ (6) 



Indeed, 6 ~ means that the particles are nonrelativistic and that the RRG is nothing but the 
dust. Furthermore, for small velocities one can just set h = (3. The main purpose of this paper 
is to establish the upper bound for the parameter h from the analysis of cosmic perturbations. 

It is easy to see that the expressions ([TJ, ([3]) and ([6]) interpolate between the dust (the 
limit 6 = 0) and radiation (6 oo) extreme cases. It is important to note that the expression 
([5]) does not represent a simple sum of the pressureless and radiation components. Instead 
this is a formula for the ideal relativistic gas of massive particles which undergoes an adiabatic 
expansion. At high red shift z — > cxd the gas is compressed and its temperature is high. In 
this case the expression above looks as ultrarelativistic one. When 2; — > — 1, the gas becomes 
almost pressureless and the above expression is close to the one for the dust. Of course, we are 
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interested in the finite time intervals and, for this reason, can not separate the RRG equation 
of state and scale dependence ([5]) to the radiation and dust parts. Due to this feature the 
expression above looks as a useful tool for various problems of cosmology. In particular, here 
we will use the RRG model as a new tool for testing the warmness of the DM today and in the 
early Universe. 

Let us consider the cosmological solution in the RRG model ^22j . The Friedmann-Lemaitre 
equation has the form 

H'{z) = ^ [p{z) + pa] + H'.Qlil + zf , (7) 

where p{z) is given by ([5]) and pA = A/SvrG. This equation can be presented in the explicit 
form 

= Hi (1 + zf + -^£L= (1 + zf ^l + b\l + zf + nl] . (8) 

Let us introduce the following two useful notations: 

^ (1 + ^)^ ff^.Pi^^ {l + z){H^)'-2nlHl{l + zf 

2,[m-^lHl{l + zY]' pt{z) [H^-nlHiil + z)^]{i-r) ' 

where the prime means derivative d/dz. In the last formulas we denoted the ratio of the square 
of the rest energy density ([2]) and the energy density p as 

r = r(z) = 4^, 
p^{z) 

and also applied the usual sum rule + + 0^ = 1. 

The cosmological model based on RRG with the presence of the cosmological constant 
admits an analytic solution [22]. This solution, as one should expect, does interpolate between 
the ones for the radiation and the dust cases. In the very early Universe, when the temperature 
is high, the evolution of the Universe is close to the one in the radiation - dominated case. At 
the other end of the energy scale, in the late Universe, the solution becomes close to the one 
for the pressureless matter case. Such interpolation between the two regimes is qualitatively 
similar to the more conventional case where the matter content is composed by a sum of the 
dust and radiation. In the conventional case, also, radiation dominates at high z and dust 
dominates at low z. However, in the RRG case we observe this property in a cosmological 
model with a single fluid. This property makes RRG model a useful tool for modeling the 
behaviour of a DM particles in the expanding Universe. The analytic zero-order solution can 
be, in principle, extended for the combination of the relativistic gas of massive particles and 
radiation, or for the combination of several distinct RRG-like fluids (for example, this can be 
done by using the method of [28\). 

3 Perturbations spectrum 

Consider the cosmic perturbations in the RRG model described above. We shall follow the 
scheme elaborated for the exploration of another model which is motivated by quantum correc- 
tions [3T] and consider simultaneous perturbations of metric, energy density and the 4-velocity 
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in the co-moving coordinates 



[/« ^ [7- + , p^p(i + s) , g^, ^ g^, + V (10) 

In the synchronous coordinates we have h^^ = 0. The perturbation of the pressure should be 
derived from the eq. ([1]), such that 



3 

In this way we arrive at the following 00-component of the Einstein equation 

, 2h h (2 - r) 

where h = dt {hkk/<^'^)- Other equations follow from the variation of the conservation law 
5 (V^Tlf) = and have the form 



6' 



4 — r / h V \ , ^ 



and 



(1 + ^) 

'P r' 5 f[\ k'^{l + z)fi 1-r 



v' + (--- + ^ S '5 = 0' 13 

where v = fi (VkSU^) and we used Fourier expansions for 5{z) and v{z) 



^(x, z) = I 5(k, z) e'"-'' , v{^, z)= J v{k, z) e**^-^ , /c = |k| 



In order to explore the equations (|TT|) . p2|) and p3j) one has to choose the initial conditions, 
related to the choice of the transfer function. We have performed the numerical analysis using 
two kinds of these functions. The more complicated one was introduced in [29] and was 
explained in details, e.g., in |30l [31]. The second option is a more simple transfer function 
from the book [7J. Both transfer functions assume a scale invariant primordial spectrum, and 
determine the spectrum today considering the Universe with the cosmological constant and 
filled by DM. Using the transfer functions we can fix the initial conditions at a redshift after 
the recombination epoch. It is remarkable that the results for the power spectrum obtained 
through these transfer functions actually coincide. The most relevant plots which show the 
comparison with the 2dFGRS observational data [2^ are presented in Fig. 1. 

The relevant quantity to be compared with the observational data is the power spectrum 
parameter defined by 

Vk = ^l , (14) 

where 5^ is the component of the Fourier transform of the density contrast S{t), which is 
computed by integrating the equations for the cosmic perturbations ([TT]) . ([T2|) and ([13|) for a 
given value of k and with a given initial conditions (as indicated above). 
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Figure 1: Power spectrum for the RRG-A model, for fixed = 0.04, f22)M ~ 0.21 and 
= 0.75 (flat Universe), with the values b = 10~^ b = W~^, b = 2x 10"^ and b = 10"^ 
The theoretic plots are presented together with the LSS data from the 2dfFGRS [27] . The 
ordinate axis represents the log of P{k) = \5m{k)\'^ at z = 0. In the abscissa we have the log of 
the wave number k given in /iMpc~^ units. 
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In the present case, since the upper bound for the possible values of b has great physical rel- 
evance, it proves useful to establish this bound with more certainty. For this end we performed 
calculation for a set of different values of b and then applied the statistical method to compare 
the result to the power spectrum data of the 2dFGRS observational program. The quality of 
the fit between the theoretical estimate and the observational data can be characterized by the 
quantity 

X^ = Ef^^^)'' (15) 



where V^. is the observational data for the power spectrum for a given wavenumber ki, V\. 
is the corresponding theoretical result obtained by the numerical integration of the equations 
for the perturbations pT|) . (fT2|) . ([13]) and cTj are the observational error bars related to the 
measurement. As smaller is the parameter x^, better is the fit. Of course, since our theoretical 
model depends on some input parameters such as 6, fi^^ and O^, the value of depends 
also on these parameters. At a first step, in order to obtain estimations for the relevant 
parameters, we reduce the three-dimensional probability distribution to the one-dimensional 
one by choosing the values of the present day cosmological parameters 

fiO^ = + = 0.04 + 0.21 = 0.25 and Vt\ = 0.75 , 

which correspond to the flat space section of the space-time manifold. 
Using the quantity the probability distribution is given by 

p = ^e-^'/2 , (16) 

where yl is a normalization constant. The final result for the one-dimensional probability 
distribution function (PDF) for the parameter b is displayed in Fig. 2. Using this plot we can 
see that this probability distribution goes sharply to zero for the values above 5 = 5 x 10~^. 
We can conclude that there is an upper bound for the "warmness" parameter 6, that means 
5 < 3 X 10-5 ~ 4 X 10"^ 

Before starting to discuss physical interpetations of our results, let us present some com- 
ments concerning the validity and restrictions of the bound on the parameter b obtained above. 
The power spectrum, which is defined in the space of the fe's, is the Fourier transformation 
of the two-point correlation function, which is defined in the real space, x. Hence, the power 
spectrum implies a Fourier decomposition. The different Fourier modes are independent only 
when the linear approximation is valid; non-linearity leads to the mixing of the different Fourier 
modes, which do not remain independent anymore. The non-independence of Fourier modes 
must be taken into account when the statistics analysis in the k space is performed. This is 
done through the use of the covariance matrix, Cij. According to this method, the likehood 
function for a given set of data Aj is given by 

£(xexp|-^AiC^iA,j, (17) 

where C~^^ is the inverse of the covariance matrix. If the conditions of linearity are satisfied, the 
covariance matrix becomes diagonal, and the likehood function reduces to the usual expression 
encoded in the parameter. 
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Figure 2: Probability distribution for the parameter b. The probabihty becomes essentially 
zero for 6 > 5 X 10~^. 

We win use the 2dFGRS data such that < fc/i"^ < 0.185 Mpc"^ [27j. This interval lies 
outside the non-linear regime, generally fixed as k h^^ ^ 0.8 Mpc^^ , corresponding to a scale 
greater than 8 Mpc. If the linear approximation is valid, the covariance matrix is diagonal, and 
the use of a statistic is justified. However, the limits on the wavelength where the linear 
regime can be safely applied is not very well established. In reference [33], it has been argued 
that the confidence on the linear approximation for the 2dFGRS data restricts the scales to 
kh^^ < 0.15 Mpc^^. On the other hand, the cosmic variance leads to a restriction such that 
kh^^ > 0.02 Mpc^^ . If we use data ranging in the interval 0.02 Mpc^^ < kh^^ < 0.15 < 
Mpc^^ the results does not change substantially, as it will be discussed latter. 

To obtain a better estimate for the velocities of the DM particles in halo of galaxies, we 
must enter deeply into the non linear regime. All considerations presented in this work are 
based on the applicability of the linear analysis and therefore we assume it to be valid. Using 
the matter power spectra data, we can not avoid in this case the use of the full covariance 
matrix, since the different modes are not independent anymore (see for example reference [34] 
for a beautiful analysis of the non-linear effects and the consequent use of the covariance matrix 
in the context of the baryonic acoustic oscillations). Hence, our previous estimations must be 
seen as a lower bound, since the process of contraction in the formation of the galaxy must 
increase the velocities of the particles. 

Furthermore, in the previous analysis, we have fixed the quantity of dark matter and dark 
energy, based on the 5-years results of Wmap. In order to verify how this restriction may 
influence the obtained bounds for the parameter b, we now consider a two-dimensional pa- 
rameter space, varying at the same time b and the dark matter quantity i^dmo (or equiva- 
lently the dark energy density 0,\o, since we are restricted to a flat spatial section). The 
results are shown in flgure 0] using the data < kh~^ < 0.185Mpc~^, and in figure [3] using 
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Figure 3: The two-dimensional probability distribution, for the flat spatial section, when both b and f^Ao are 
varied (left), considering the modes < kh~^ < 0. 185 M pc~ ^ . Higher probabilities are indicated by the brighter 
regions. The one-dimensional probability distribution for b, after marginalizing on Qao is shown at the plot on 
the right. 

0.2Mpc~^ < kh'^^^ < 0.150Mpc~^. According to these plots, the bounds for the value of b 
remain essentially the same when we perform these additional variations of other cosmic pa- 
rameters. Moreover, they are also in a good agreement with the case when we used only the 
best fit value for Q\q and ^dmo given by Wmap. 

Now we are in a position to discuss physical significance of the bounds obtained from 
cosmic perturbations analysis. The restriction for the parameter b which we derived from the 
numerical analysis of cosmic perturbations, can be easily translated into the bound for the 
average velocity (or warmness) of the DM constituents. For this purpose we have to note 
that b is necessarily small and therefore the relation ([6]) converts into b = (3 = v / c. Then we 
arrive at the bound for the average speed of the massive relativistic particles of the DM in 
the present-day Universe v < vq, ~ 10 — 12 km/ s. This bound agrees with the standard 
evaluations (see, e.g., [6j 4]) obtained from the numerical simulations of the structure formation 
in the neutrono-dominated Universe and also from the model-dependent considerations. 



It looks tentative to compare this bound for the speeds of the DM particles with some 
astronomic observable, as it is usually done in the framework of standard approaches [26]. One 
can, for instance, try to compare this bound to the known one for the spiral galaxies, which is 
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Figure 4: The two-dimensional probability distribution, for the flat spatial section, when both b and f^Ao are 
varied (left), but restricting the modes to 0.2Mpc~^ < kh~^^ < 0.150Mpc~^. Higher probabilities are indicated 
by the brigther regions. The one-dimensional probability distribution for b, after marginalizing on Q.ao is shown 
at right. 

aboulEI 2A0 km/ s. Obviously, there is no correspondence between the two numbers. However, 
let us note that the galaxy is an object which definitely lies out of the linear perturbation 
regime which we deal with here. In general, any comparison of the results obtained from the 
linear approximation to the cosmological perturbations, including the one for the upper limit of 
the velocity of DM particles with the dynamics of virialized systems in the Universe (galaxies, 
clusters of galaxies, etc.) must be performed with great caution, because those virialized 
objects are in the deep non-linear regime, with typical densities hundreds of times larger than 
the critical density. 

In general, the estimations of velocity of dark matter particles strongly depend on the 
nature of the DM candidate under consideration and especially on the reference frame where 
the velocity is evaluated. For example, in the paper [35] were obtained a dispersion of primordial 
velocities (which is not affected by non-linear effects) with respect to the Hubble flow. The 
result is of the order of magnitude about 10^^^ in unities of the velocity of light for the 
WIMPS with a typical mass of a few GeV^s and about 10^^'' in the same units for the axion, 
even if the typical mass of the axions is some 10^^ orders smaller than the typical mass of the 
WIMPs. This case shows how the nature of the dark matter candidate influences the estimate 
for its primordial velocity. In reference [36], the possibility that the warm dark matter can be 
described by sterile neutrinos has been analysed. The possible mass range goes from 1 keV 
to tenths of keV. Now, using the results for the WIMPS quoted above, and the usual non- 

^This approach is close to the one of [32], which is based on the Maxwellian dark matter velocity distribution 
for spherically symmetric and isotropic halo. 
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relavistic relation for the ratio of velocities of particles having the same energy but different 
masses, we find an upper limit for the velocity of the warm dark matter particle of the order 
of 10~^ in unities of c. This is far below our upper limit near the maximum of probability 
distribution. On the other hand, the two estimates correspond to distinct reference frames and 
therefore these results do not contradict our analysis. 

We could try to use, for instance, the data of weak lensing investigations concerning the 
dark halo of galaxies. But, the results obtained in this way still concern a very non-linear 
regime, see for example |37j . In order to chose the appropriate object, we note that among the 
galaxies, the dwarf spheroidal elliptical ones are those with the large proportion of dark matter 
to baryonic matter, with a mass/luminosity ratio that can be as large as 500, in solar units, 
while ordinary galaxies have a ratio of the order of some 10 — 50 [35]. Even if these objects are 
deep in the non-linear regime it is remarkable that the dispersion velocity becomes essentially 
constant far from the center with a typical value of the order of v ~ 10 km /s [39], comparable 
with the velocity bounds we have obtained for the dark matter particle. 

Since the dwarf spheroidal galaxies represent the extreme case of virialized system domi- 
nated by dark matter, such agreement of the typical velocities with our bound is perhaps not 
meaningless, even if a clear determination of the extension of the dark halo would be necessary 
to put this comparison into more solid grounds. For the dwarf irregular or spiral galaxies the 
typical radial speed of the stars is evaluated to be about 10 — 12 km /s [40j and we meet a 
nice correspondence with our bound. It is amusing that we arrived at this correspondence by 
using a very simple RRG model [22] and not the complicated approach based on the Einstein- 
Boltzmann system [26j . 

The last problem which we can easily address within the RRG model is the dynamics of 
the DM average speeds in the expanding Universe. In other words, it would be interesting to 
calculate how this speed depends on the red-shift parameter z, or on the temperature Tcmb 
of the cosmic background radiation (CMB). In order to address this issue one has to use the 
relations ([6]) and ([5]). The unique role of the parameter b is to define the "warmness" of the 
DM in the last of these relations, so it is obvious that at the higher z we have h{z) = 6(1 -|- z), 
where h is the modern value. Furthermore, since the upper bound for h nowadays is much less 
than 10"'^, for the potentially relevant z < 1000 we can, according to ([6]), safely use the formula 
V = be. In this way we arrive at the following relation for the average speed of the DM particles 

v(z) = cb(z) = cb(l + z) = V X — , (18) 

^CMB 

where we used Tcmb ~ (1 + -2) and denoted T^lj^, b and v the corresponding quantities 
for z = 0. So, in the framework of the RRG model we note that the average kinetic energy of 
the DM particles is proportional to the square of the CMB temperature, the result which is 
familiar from the conventional (but more complicated) considerations (see, e.g., [9J)- 

4 Discussions and conclusions 

We have considered the structure formation in the model where the DM is described by the ideal 
relativistic gas of identical massive particles. Instead of using Maxwell distribution, we have 
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employed the RRG model [22] which is closely reproducing Maxwell distribution and, at the 
same time, is rather simple. As a result we arrive at the strong limit on the parameter b, which 
should satisfy the upper bound b < 3— 4x 10~^. According to the relation ([6]), this is equivalent 
to the upper bound on the velocities of the DM particles v < vq = 3 — 4 x 10~^c = 10 — 12 km/ s. 
This restriction is much more severe than, e.g., the one discussed earlier in |4 H 142 1 1 ^ I44j on 
the basis of the nonrelativistic Maxwell distribution [42j and is essentially smaller than the 
typical velocities of the stars in the spiral galaxies. Also, it is about two order of magnitude 
smaller than the escape velocities for the spiral galaxies. 

Does our result mean that the actual velocities of DM particles can not be greater that the 
mentioned bound vqI An obvious answer is no. Let us remember that both DM and baryonic 
matter can acquire an extra kinetic energy after the galaxy starts to form and the linear regime 
of the cosmological perturbations can not be applied. One can see the corresponding process 
as a kind of the usual transformation of the potential gravitational energy into the kinetic one. 
This process has nothing to do with the linear perturbations we have studies here. However 
taking smaller astrophysical objects such as dwarf galaxies, we arrive at the surprisingly nice 
correspondence between the observed average speeds of the stars in such galaxies and our upper 
bound vq. This correspondence shows that the RRG model is, perhaps, the simplest way to 
arrive at the reasonable estimates concerning not only the behavior of the conformal factor, 
but also the linear cosmological perturbations. 

The result described above is universal in the sense there is no dependence on the origin 
and properties of the WDM constituents. In particular, the restriction on velocities does not 
interfere with the one for the masses of the WDM particles, which can be even macroscopic ones. 
We were just treating them at a component of ideal relativistic gas and derive restrictions on 
their velocities. Hence, these restrictions apply equally to all known models of WDM, see, e.g., 
refs. [45 1 [2H [T8| [20l I46j . For example, they apply to the models of DM particles which do not 
interact with other matter and with themselves, except gravitationally [8| 1471 According 
to our results, even this simplified model can produce good predictions for the spectrum of 
linear perturbations, if we assume the ideal gas equation of state for these weakly interacted 
particles. Another interesting point is that our results show that the anisotropy of the DM 
distribution does not play a critical role in the definition of the perturbations spectrum and 
that the equilibrium distributions such as the Maxwell one (which is closely reproduced by 
RRG) is, in principle, sufficient to arrive at the reasonable estimates for the speed of the DM 
particles. 

Last, but not least. The RRG model may be successfully applied for the investigation of 
more complicated situations, including two distinct non-interacting ideal gases |49j and also 
may be useful in describing the interaction between these gases, e.g. through the use of viscosity 
(see, e.g., [25]). We expect to explore these issues in the near future. In general, our model 
proved useful in exploring relativistic properties of the ideal gas of massive particles, it can be 
applied for solving various problems of gravitational physics. 

In order to verify to which extent the model described above may correct the excess of 
power of the ACDM at small scales, the non-linear regime may be explored. As a result one 
may hope to achieve a more detailed description of the structure formation. At the moment it 
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is unclear whether the RRG model can be a useful tool in such case. We hope to explore this 
issue at the consequent stage of our work. 
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